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The equations of motion of a gravitationally interacting two point-mass system have been derived some years ago 
/"■n \ up to the 5/2 post-Newtonian (2.5PN) approximation^], in harmonic coordinates Recently, it has been possible 

■ to derive the third post-Newtonian (3PN) Hamiltonian of a two point-mass system [HI within the canonical formalism 
of Arnowitt, Deser and Misner (ADM) ||. It was found that, at the 3PN level, the use of Dirac-delta-function 
sources to model the two-body system causes the appearance of badly divergent integrals which (contrary to what 
happened at the 2.5PN ||,|| and 3.5PN j?]] levels) cannot be unambiguously regularized [§|,|],[§. The ambiguities in 
the regularization of the 3PN divergent integrals are parametrized by two quantities: w s tatic and Wkinetic- 

Prompted by a recent remark fiofl , the purpose of this work is to show that requiring the (global) Poincare invariance 
[ of the 3PN ADM Hamiltonian dynamics uniquely determines one (and only one) of these regularization ambiguities: 
• namely the "kinetic ambiguity" parameter akinetic- [The "static ambiguity" w s t a tic remains unconstrained because 
it parametrizes a 0(c~ 6 ) Galileo-invariant additional contribution to the 3PN Hamiltonian.] Parallel work in the 
harmonic-coordinates approach to 3PN dynamics has recently obtained similar results . 

Note that general relativity admits (when considering isolated systems) the full Poincare group as a global symmetry. 
Therefore, whatever be the coordinate system used (as long as it respects asymptotic flatness) the general relativistic 
^ dynamics of A^-body systems should embody some representation of this global Poincare symmetry. When solving 

y| Einstein's equation by a weak- field, "post-Minkowskian" expansion, ■ s fgg iLV — r\^ v = = Gh^ + G 2 /i™ H , and 

fc"p." fixing the gauge by the "harmonicity condition" , d v = 0, the whole scheme stays manifestly invariant under the 
usual (linear) representation of the Poincare group: x'^ — A£J x v 4- a M (assuming that the regularization procedure 



used to deal with the point-mass divergencies is manifestly Poincare invariant). In such a case the A-body dynamics 
will be invariant under the representation of the Poincare group induced on the dynamical variables, say x l a (t), 
x l a (t), a = 1, . . . , A, by the action of the usual linear Poincare transformations. This global Poincare symmetry has 
been explicitly checked at the 2PN level in Ref. jlj] by proving that the 2PN (acceleration-dependent) two-body 
Lagrangian in harmonic coordinates 1 13 |] changed only by a total time derivative under a generic, infinitesimal 



Poincare transformation. In this work we consider the 3PN two-body Hamiltonian derived by Ref. |4| within the 
ADM canonical formalism. This formalism is not manifestly Poincare invariant because it splits space and time, and 
fixes the coordinates by the following gauge conditions: S^tt^ = 0, dj (cjij — ^g ss 8ij) = 0. This lack of manifest 
Poincare invariance is not problematic (though it introduces some technical complications). Indeed, we shall explicitly 
show in this paper that the global Poincare symmetry of the two-body dynamics can be realized in phase space, albeit 
by a somewhat complicated, nonlinear action. 

The basic principle that we shall follow to study Poincare invariance of the 3PN two-body Hamiltonian if(x Q ,p a ), 
a = 1, 2, with its associated Poisson bracket structure, 



{ A(x Q ,p a ),s ( x a ,p Q)} (D 



dA dB dA dB 



1 We recall that the "nPN approximation" refers to the terms of order (v/c) 2n ~ [GM/(c 2 r)) n in the equations of motion. 
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is the following: the presence of a Poincare symmetry is equivalent to requiring the existence of "generators" P M , J^ v , 
realized as functions P^(x a ,p a ), J M1/ (x a ,p a ) on the two-body phase-space (xi, X2, Pi, P2), whose Poisson brackets 
(||) satisfy the usual Poincare algebra (here we set c = 1): 

{P^,P'} = 0, (2a) 
{P 11 , J pa } = -^ p P a + ^ P p , (2b) 

{ J"", J pa } = -rf" J pa + rf p J va + if* J pv ~ rf v J PA \ (2c) 

where r]^ = diag (-1, +1, +1, +1). 

The functions P M (x a ,p a ), J Miy (x a ,p ) generate (in phase space) the infinitesimal Poincare transformations: 
Sa,u F = {F, a p P^ + \ tji p ' v J^}. Finite transformations are then (in principle) defined by exponentiating these 
infinitesimal actions. The satisfaction of the algebra (||) ensures that one thereby generates a consistent Poincare 
symmetry. The time component P° (i.e. the total energy) is realized as the Hamiltonian P(x Q ,p a ) (including the 
rest-mass contribution). The other generators can be decomposed as: P l (three momentum), J % = i e lkl Jke (angular 
momentum), and K l = J t0 (boost vector). One further decomposes the boost vector K l (which represents the con- 
stant of motion associated to the center of mass theorem) as K % (x a , p a ; t) = G % (x a , p a ) — t P l (x a , p a ) so that the total 
time derivative dK l /dt = dK z /dt + {K 1 , H} = —P l + {G 1 , H} = 0. Finally, the Poincare algebra explicitly reads 

{Pi,H} = {Ji,H}=Q, (3a) 

{ Ji, Pj} = Eijk Pk , {Ji, Jj} — £ijk Jk , (3b) 

{Ji,Gj} — EijkGk 7 (3c) 

{Gi,H}=Pi, (3d) 

{G i ,P j } = ±H6 ij , (3e) 

{Gi,Gj} — —— e^k Jk ■ (3f) 

As the gauge fixing used in the ADM formalism manifestly respects the Euclidean group (which means that 
P(x a ,p a ) is translationally and rotationally invariant), the generators Pi and Ji are simply realized as 

P(x Q ,p a ) = ^Pai , (4a) 

a 

Ji(x Q ,p a ) = ^2e iki XaPae, (4b) 



and exactly satisfy Eqs. (3a) and (pq). The condition fl3c| ) will also be exactly satisfied if Gi is constructed as a 
three- vector from x and p a . Finally, the condition for full Poincare invariance boils down to the existence of a vector 
Gi(x a j_p a ) satisfying the three non-trivial relations (|3c|), ( |3c| ) and (pq), in which enters, besides Pi and Ji given in 
Eqs. (|4a| ) and (|4q ), the full (3PN- accurate) Hamiltonian: 

H(x a ,p a ) = ^2m a C 2 +P N (x a ,p a ) + i-Pi PN (x a ,p a ) + -ji72PN(x a ,Pa) + -g #3Pn(x , Pa) + O (-) . (5) 

At the Newtonian order, i.e. when keeping the rest-mass term S a m a c 2 and the Newtonian-level Hamiltonian, 

P N (x a ,p a )=£— , (6) 

(r ao = |x a — X(,|), it is easily checked that the usual Newtonian center-of-mass vector 

G]v(x a ,p a ) = ^m a x l a (7) 
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satisfies Eqs. (3d)-(|3|). [Note that, in this approximation, the right-hand side of Eq. (|3eJ) yields {J2 a m a) from the 
rest-mass contribution to H.] 

To study the existence of G % beyond the Newtonian approximation, we need the explicit expressions of the 1PN, 
2PN and 3PN contributions to the Hamiltonian (ph in an arbitrary reference frame. The 1PN contribution, 
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has been known for a long time. The operation "+(l < — > 2)" in Eq. (^) denotes the addition for each term in Eq. (||) 
(including the ones which are symmetric under label exchange) of another term obtained by the label permutation 
1 < — ► 2. The 2PN-accurate explicit expression of H (x a , p a ), in the ADM formalism, was derived in Ref. flj (Eq. (2.5) 
there). [The corresponding explicit Lagrangian L^p^(x a ,x a ) is given in Ref. jL5|.] These results corrected earlier 
results by Ohta et al. p6|. The final result reads: 



i? 2 PN(x a , Pa) 
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Ref. Q derived the 3PN-accurate ADM Hamiltonian restricted to the center- of -mass reference frame: Pi + P2 = 0. 
For the present work we have generalized Ref. Q in deriving 7?3pn in an arbitrary reference frame. Our starting point 
for doing this calculation is the improved form of the 3PN Hamiltonian, i?3PN, given in Appendix A of Ref. § (Eqs. 
(A8)-(A10) there). Note first that i?3PN defined there denotes the higher-order Hamiltonian iJ3PN(x a , p a , x a , p a ) de- 



fined by eliminating the field variables h[j , h[j in the "Routh functional" i?(x a , p Q , , ) introduced in Eq. (33) 

of Ref. B. However, it was shown in Ref. Q that one could reduce the higher-order Hamiltonian _ff 3 pN(x a , p a , x a , p a ) 
to an ordinary Hamiltonian iJ3PN(x^, p^), at the price of the following (3PN-level) shift of phase-space coordinates: 



dH 

dpa 



Pa = Pa 



dH 
<9x Q 



(10) 



After performing the shift (|10|) with respect to the original ADM coordinates x a , p a (we henceforth drop the primes 
for notational simplicity), the calculation of the 3PN (order- reduced) Hamiltonian consists in evaluating three very 
complicated integrals: 



#3PN = —77^ 



128 ^ 
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d A x(h\ 



red 



(ii) 



The integrands hi, hi, h% are given in Eqs. (A9) of |J. The order-reduced integrand /i^ cd is defined (as shown in 
by using the Newtonian equations of motion to eliminate x a and p a when computing the time derivative hf^ (which 
enters the last two terms of hi). As explained in g], this new form of the 3PN Hamiltonian is free of "contact term" 
ambiguities, and the integrals it contains can all be uniquely defined by using the Riesz-type regularization procedure 
explained in [Q. We have recomputed from scratch all the integrals by using the generalized Riesz formula given in 
[Q| . This Riesz-regularized 3PN Hamiltonian reads explicitly (in an arbitrary reference frame )0 



The authors thank Luc Blanchet for pointing out that two terms were missing in the printed version of our non-center-of-mass 
Hamiltonian. 



3 



^3PN^ a ,P„j- ^ — 



..(Pi) 3 , , ((Pi-P2) 2 + 4p?p!H , (Pfpi-2(P1- P2 ) 2 )(pi- P2 ) 



1Q (Pi ("12 • P2) 2 + pj (ni2 • Pi) 2 )p 2 24 P 2 (Pi -p 2 )(ni2 -Pi)(ni 2 ■ p 2 ) g pf (p i • p 2 )(ni 2 • p 2 ) 2 



4 2 



4 2 



3 3 



(7p 2 p| - 10 (pi •p 2 ) 2 )(ni 2 -pi)(ni2 -p 2 ) p 2 (n 1 2 -pi) 2 (n 12 • p 2 ) 2 



3 3 



4 2 

1 'JTl 



+15 (Pi • P 2 )(n i2 ■ Pi) 2 (ni2 ■ p 2 ) 2 _ 18 P 2 (ni2 • Pi)Qii 2 ■ p 2 ) 3 + g (nti 2 • Pi) 3 (ni 2 • p 2 ) 3 



to 3 to 2 



to 3 to 2 



m\m\ 



G 2 mim2 



1, 07 > 2 ) 2 115 P 2 (Pi-p 2 ) , 1 25( Pl -p 2 ) 2 + 371p 2 p 2 
(mi - 27m 2 )— 34 — TO i a 1" — m 2 " 



16 



16 mfm 2 48 



17p 2 (ni2 ■ pi) 2 1 (l5p 2 (ni2 • p 2 ) + 11 (pi • p 2 ) (ni 2 • pi))(ni2 • pi) 5(ni 2 -pi) 4 
■ mi- 



16 mf 



m 3 m 2 



12 mf 



3 (ni2 • pi) 3 (ni2 • p 2 ) 125 (pi • p 2 ) (n i2 • Pi)(ni 2 • p 2 ) 10 (n i2 • pi) 2 (ni 2 • p 2 ) 2 
-mi § h -rrr^ h — m 2 

2 mjm 2 12 tojto^ 3 m^m. 



1 Aion i in-3 \Pl( n 12 • P2) 
- — (220mi + 193to 2 ) 3 — 5 

48 mfrnj 



G 3 mim 2 



— (^466 to? + (473- ^7T 2 Wm 2 + lSOrniA^ 
48 \ V 4 / / to 2 



( 77(m 2 + to 2 ) + f 143 - ^7t 2 )toiTO 2 ) ^ Pl ' P2 ^ + (61 to 2 - (43 + ~y ) /// , m • 
16 \ \ 4 / J TO1TO2 16 \ V 4 

1 /„/ 2 2s , fna , 3 2 \ \ (1112 -pi)(ni2 -p 2 )] 

( 21(to>i + to 2 ) + |119 + —it jm 1 m 2 J- 



\ (ni2 • Pi) 



toito-2 



1 G 4 toito1 



227 21 

— — 7T j /;/ I -!■ 1)1-2 



(1^2). 



(12) 



However, it was emphasized in Refs. ||,||,|| that the nature of the divergent integrals which had to be regularized 
to compute H^p N , Eq. (|l2|), was such that the result should be considered as being partly ambiguous. These regu- 
larization ambiguities have been discussed in Refs. 0], and, in more detail, in the Appendix A of B. We have 
recomputed, in an arbitrary reference frame, the various regularized versions of all the momentum-dependent formal 
"exact divergences" A31, A32, . . . , A38 defined in Appendix A of ||. These contributions should formally vanish, 
but their regularized values do not vanish and thereby exhibit the regularization ambiguities present at 3 PN. We 
find (in confirmation of the result given in the Introduction section of Ref. B) that all the momentum-dependent 
regularization ambiguities are equivalent to adding to Eq. ( |l2| ) a term of the (specific^) form 



#3 k ra ic (xa, Pa) = +1 kinetic ^ — [p 2 - 3(n 12 • Pl ) 2 + p 2 - 3(n ia • p 2 ) 2 ] 



(13) 
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where ^kinetic is an arbitrary parameter. In addition to this "kinetic" regularization ambiguity it was pointed out in 
H and H, that there is also a "static" (i.e. momentum-independent) regularization ambiguity of the form 



HlvN C ( x a, Pa) 
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Note in particular the absence of terms mixing pi and P2. 
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where w sta tic is a second arbitrary parameter. Finally, the 3PN (order-reduced) Hamiltonian is of the form 

-ff3PN(x a , Pa) = #3PN + #3PN tlC + #3PN C > 

and depends on two, up to now undetermined, real parameters osmetic and astatic- 

The problem to solve is now the following: does there exist a (3PN-accurate) center-of-mass vector, of the generic 
form, 

G 4 (x a , p Q ) = ^(M a (x k , p,,) < + A a (x b , p 6 ) pi) , (16) 

a 

where M a and N a are scalars that reduce to m a and 0, respectively, in the Newtonian approximation, such that 
Eqs. (|3dl )-(j3fl) are fulfilled (within the 3PN accuracy) when the Hamiltonian is given by inserting Eqs. (||), (|J), (j|) 
and ( p_5| ) in Eq. (||)? We have tackled this problem by the method of undetermined coefficients, i.e. by writing the 
most general expressions for the successive PN approximations to the functions M a (xb,p ) and iV a (x6,Pf>), 

M a =m a + -^M a +-^ M a +^ M a > N a = -^N a + -g N a , (17) 

as sums of scalar monomials of the form: c„ „ 1 „ 2 „3„ 4 „ 5 rj" 2 ™ (p?)™ 1 (p^)™ 2 (pi ■ p 2 ) n3 (ni 2 ■ Pi)™ 4 (n i2 • p 2 )™ 5 , with 
positive integers rio, . . . , rig. Besides dimensional analysis (which constrains the possible values of no, . . . , ng at each 
given PN order), and Euclidean covariance, including Parity symmetry, wc only required time reversal symmetry 
(which imposes that M a be even, and N a odd, under p a — ► — p a ). We did not impose any a priori constraints on the 
mass dependence of the coefficients c n (m\, m^), nor did we use the 1 <-> 2 relabeling symmetry. 
The 1PN approximation to G % being well-known (see, e.g., p7||), 

m ipn = 1p! _lgmim a 
2 mi 2 ri2 

jVf PN = 0, (18b) 
with M2 PN obtained by a 1 < — > 2 relabeling, we started looking for the most general G 1 at the 2PN level. At this 



level, there are 20 unknown coefficients c n , and Eq. (3d) yields 40 equations to be satisfied. We found that there is a 



unique solutionn to these redundant equations, namely 
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with M| PN and ATf PN obtained by a 1 < — > 2 relabeling. 

We have a posteriori checked that this unique ADM-gauge, 2PN ccnter-of-mass vector agrees (after taking into 
account the shift x^ DM = z a — S*z a (z, z) [fL4| ) both with the harmonic-gauge 2PN G l (z a ,z a ) first derived in Ref. 
p2[ , and with the Landau-Lifshitz-likc |Q, ADM-gauge calculation of G*(x a ,x a ) performed in Ref. ]l8| ]. We have 
also checked that the remaining Poincare-symmetry constraints, Eqs. ([k]) and (Q), are also fulfilled. Concerning Eq. 
(pa), it is easy to see, in general, that it is equivalent to the constraint 



E 

a 



M a (x b ,p b ) = ^H(x b ,p b ). (20) 



At the 3PN level, the most general ansatz for M^ PN , N^ FN , involves 78 unknown coefficients c„, while Eq. ( |3d| ) 
yields 138 equations to be satisfied. The quantity ^kinetic parametrizing the momentum-dependent regularization 



All the algebraic manipulations reported in this paper were done with the aid of MATHEMATICA. 
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ambiguity (|13|) in the 3PN Hamiltonian enters the system of equations for the unknown c„'s. [Indeed, it was recently 
noticed that H^p^ tlc is not separately boost-invariant flltf -l By contrast, the other regularization ambiguity ([l4|) drops 
out of the problem (because H|pn 1c is Galileo invariant). We found that there was a unique value of akinetic for which 
the system of equations to be satisfied was compatible, namely 

41 
24' 



^kinetic 



If ^kinetic 7^ 41/24, the 3PN Hamiltonian does not admit a global Poincare invariance. If akinetic 
unique solution to Eq. (|3d|), namely 
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We have then checked that this unique solution does satisfy the remaining Poincare-symmetry constraints, Eqs. (pa) 
and (^|), or, equivalently, Eqs. ( p(i| ) and (H). It is to be noted that the last two momentum-dependent terms in Mf , 
proportional to (1112 • pi) 2 /mi — (1112 • p2j 2 /m2, are antisymmetric in the labels 1 < — ► 2 and therefore drop out in 
the constraint (p0|), which reads Mf PN + M| PN = i?2PN- In fact, the corresponding monomials appear nowhere in 
H2PN, but must crucially be included in M^ PN . 

The main conclusion of this work is therefore that the necessary existence of a global Poincare symmetry in the two- 
body problem uniquely fixes the regularization ambiguity parameter ^kinetic to the value ( pl| ) . The explicit realization 
of this Poincare invariance is then defined by the phase-space generator G l (x a ,p a ) defined by Eqs. (|l6|), (|l7|), Jlij|), 
©,and@. 

Within the ADM formalism it would be very difficult to implement a Poincare-invariant regularization procedure. 
(The situation is different in harmonic coordinates, where one can conceive a Lorentz- invariant regularization Jil"|.) 
It is very satisfying (and in keeping with the general lore about renormalization theory) that we were able to use a 
non-Poincare-invariant regularization, but then, a posteriori, correct for it in a unique way. There remains, however, 
a last regularization ambiguity^, Eq. (fTi]), which has all the needed global symmetries and cannot be fixed in this 
way. 
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